More on scattering of Chern-Simons vortices 

Jacek Dziarmaga 
Jagellonian University, Institute of Physics, 
Reymonta 4, 30-059 Krakow, Poland * 

December 20, 1994; revised January 13, 1995 



in 
m 



Abstract 

I derive a general formalism for finding kinetic terms of the effective Lagrangian for slowly moving Chern-Simons vortices. 
Deformations of fields linear in velocities are taken into account. From the equations they must satisfy I extract the kinetic 
term in the limit of coincident vortices. For vortices passing one over the other there is locally the right-angle scattering. The 
method is based on analysis of field equations instead of action functional so it may be useful also for nonvariational equations in 
nonrelativistic models of Condensed Matter Physics. 
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I p ^ Introduction 

TVTTijr p is quite a large family of field-theoretical models both relativistic |[ |[ ||, U, ^, |l^, ^ and nonrelativistic j^, ||, |l^, ^ with 
soli|ons which posses Bogomol'nyi limit. In this limit static field equations in a given topological sector can be reduced for a 
mmimal energy configuration to first order differential equations. These equations generically admit static multisoliton solutions 
'c^c^acterised by a finite set of parameters such as positions of solitons and their internal orientations. The configurations are static 
- ifave think about solitons as about particles there are no net static forces between them. 

it was an idea of Manton that low energy scattering of monopoles in the Bogomol'nyi-Prasad-Somerfield model ^ can 
be modeled by reduction of the dynamics to a finite-dimensional manifold of parameters of static multisoliton solutions. The kinetic 
part of the Lagrangian of the original theory after integrating out spacial dependence of the fields with by now time-dependent 
parameters yields the kinetic part of effective Lagrangian quadratic in time derivatives of parameters. A metric on the moduli space 
can be read out of it. A fundamental idea in this approach is that configurations satisfying Bogomol'nyi lower bound are at a bottom 
of a potential well. A slow motion of solitons can lead only to small deformations of fields with respect to static configurations. 



The idea was succesfuUy appfied to scattering of monopoles, vortices in Abelian Higgs model 17, |T6|, CP^ solitons 

1 18, 1. Recently also extensions of the method to the case of Chern-Simons vortices both relativistic and nonrelativistic were done 

1 19, 2^, ^ |2^. However as was first pointed out in in these cases a new problem arises. Lagrangians of these models 
contain terms linear in time derivatives such as Chern-Simons term and/or Schrodinger action. By just promoting parameters to 
the role of collective coordinates one can reliably calculate only terms in the effective Lagrangian linear in velocities. To compute 
kinetic part one has to take into account small deformations of the fields with respect to static configurations. It is enough to 
consider only deformations linear in velocities. Such deformations can be in principle calculated from full field equations linearised 
in deformations and terms linear in velocities. However also terms linear in accelerations and third time derivatives arise and as I 
have discussed in there is no appearent reason why they could be negligible as compared to velocities. Such an "approximation" 
can lead to serious inconsistencies. 

In this paper I put the problem on a slightly different footing. The acceleration terms are not neglected. There are no net static 
forces between vortices so their accelerations must be zero for vanishing velocities. Thus we can assume that acceleration vector is 
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at least linear in velocities. We can expect this linear term to be nonzero because there are charge-flux interactions by Lorenz-like 
forces. Thus acceleration is not neglected but replaced by a position-dependent matrix lo times velocity. The same procedure can 
be applied iteratively to the third time derivative. Finally in a special limit of coinciding vortices such a unique form of cj-matrix 
is extracted which admits regular deformations of fields. A knowledge of uo is enough to establish the form of moduli space metrics. 
A brief comment on how this approach works in Abelian Higgs model is added. 



2 The model, zero modes and useful notations 

We take the Lagrangian of the self-dual Chern-Simons-Higgs model in the form 



L = + D^(j)*D^q^ - V{\ ct> I) , (1) 



where V{\ |) = ^ | |^ (| P — D^(f) — d^cf) — A^(f>, the signature of the flat 2+1 dimentional metrics is — , — ) and the 
Levi-Civita symbol is choosen so that e"^^ = — 1. A variation of the action with respect to Aq leads to Gauss' law constraint 

20Xatx - ^o) = , (2) 

where we have introduced (f) ~\ cj) \ e*^ and the magnetic field is i? = — -Fi2- When one takes into account the Gauss' law one can 
effectively rewrite the original Lagrangian in a new form 

We have just introduced f —\ (f>\ and such that e^^ = 1. The energy density for a static configuration with a positive topological 
index is 
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e = [9,/ - fe.k{dkX - Ak)Y + + — —Y -B . (4) 



The magnetic flux is quantised as 27rn with n being the winding number. Thus energy is bounded from below by this value of 
magnetic flux. If k > this Bogomol'nyi lower bound is saturated by static configurations being solutions to first order differential 
equations 

dif = fSikidkX - Ak) , (5) 
s,,d,A, = ^p{l-p) , (6) 

1 - 

^0 = ■ (7) 

These equations admit static multivortex solutions parametrised by a set of 2n real parameters. In Coulomb gauge one can take 
X = Arg{z — Zy), where the sum runs over vortices labeled by w's and complex coordinates on the plane are used. In this 

gauge there is only one second order differential equation to solve 

V2ln/ = 4/2(/2-l) + 2^^,5(2)(z-z,„) . (8) 

v — l 

Once / is known other fields can be obtained from Eqs.(||j|). The singular sources on the R.H.S. of the above equation enforce 
modulus / in a close vicinity of p-fold zero zq to behave like | z — zq |p (it is a leading term of the expansion). This equation in 
particular admits cilindrically symmetric vortex solution with winding number n (p ^ f[r) e-xjpin9. We will parametrise multivortex 
solutions in following equivalent forms 

n-l 

4>={z- zi)...{z - Zn)W{z, z, z,) = (z" - ^ \kz'')W{z, z, z„) (9) 

fc=0 

and 



iinfffi;^-^-) =ii,^"-i:'£o/^^ 



x = ^in ^'r'): 7\ ^^^ .. r^^iT., ■ (10) 



A's are complex coefficients of the n-th degree polynomial with roots z^, and is a positive real function. Although the parametri- 
sations are equivalent it will later appear that A's are more efficient in description of vortices passing one over another. 
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Each of these multivortex solutions posesses 2n zero modes. 

SAk = dkSx + ekidi^ ^{dk-^+ekidih^)SX^ , 

Mo = -^fhfSX^ , (11) 

Hi 

where Ap = + iX^ and — (1, i). Equation (||) hnearised in fluctuations becomes 

\/^h^ + - 2p)/^^ = . (12) 

p denotes moduli squared From the asymptotics of / close to its zeros we can extract the leading term in fluctuations 

~ -Re{—^^ ) (13) 

as (z" — Xgz'^) ~ 0. This is the only singular term in the expansion around the actual zero of the Higgs field. This singularity is 
fine-tuned by the singularity in 6x to yield regular SA^ (see (|lT|)). 

Now because of future applications let us take a closer look at the coincident n- vortex solution (p = f{r)eKpin9 with /(r) 
satisfying 

If" + — -{ff = \l\f~l) (14) 
r 

with boundary conditions /(O) = and /(oo) = 1. Close to zero it behaves like / ^ /or" — 2K^{n+iy r^"^"^ + ■■■■ Fluctuations for 
small A's can be written as 

h{r, 9) = hf{r, 0)X^ = -Hn^p[Xl,_p cos(n - p)0 + A^^ sin(n - p)9] . (15) 

H's satisfy following equations 

^T2+-T' ^-^)Hn-p + -P(l - 2p)i/„_p = , (16) 
dr^ r ar 

with a normalisation that close to zero Hm ^ -k^ and it asymptotically vanishes at infinity. 



3 Slowly moving vortices 

For more clarity in this paragraph and in what follows we will rescale gauge fields A^ — > k~^A^ and coordinates x'^ k~^xP'. On 
the level of field equations it amounts to fixing k ~ +1. 

The aim of this paper is to investigate slow motion of vortices in adiabatic approximation. In the case of Chern-Simons solitons 
a new difficulty arises because there are terms linear in time derivatives in the Lagrangian of the theory. Only terms linear in 
velocities can be correctly calculated by direct application of former methods. As was pointed out in |20| to obtain the kinetic term 
one has not only to make the static fields time-dependent by time-variation of their 2n parameters but one also has to take into 
account deformations of the "static" fields up to terms linear in velocities. It is something like a generalised Lorenz transformation. 
In this paragraph we will develope general formalism for calculating such corrections. It will appear that to have regular solutions 
one can not neglect terms proportional to accelerations as it was anticipated in [ p2[ . 

Let us make the fields time dependent by variation of parameters and add to them deformations linear in velocities, say 
/(z, A) f[z, X{t)] + Af[z, X{t),X{t)] etc. Evaluation of the effective Lagrangian goes by substitution of such fields into Lagrangian 
(^) and integrating out their planar dependence. The only terms which can contribute to the part of the effective Lagrangian linear 
in velocities are 

L(y^= J <fx[Bx+\e'=A,A,] . (17) 

and one needs to take into account only effects of promoting parameters to the role of collective coordinates. All the other terms 
are quadratic in velocities. The contribution of the second term in Eq.(p7|) vanishes in Coulomb gauge. For fairly separated vortices 
one can approximate under the integral B = — 27r^^ 5^'^^{z — z^) to obtain 

L^^ « -2^^ ^ Arg(z„ - z^) . (18) 
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With the exphcite form of the phase ( [l0| ) and some integration by parts one can rewrite the hnear Lagrangian as 

4// = -2vr^i;A,(z,) , (19) 

V 

where Zy = xl + ix^ . 

The second order term of the effective Lagrangian is 

4// = /' + 2/A/ + A/A/-(a,A/)2-ly"(/)(A/f 

+ fix - AAo)2 - f{AA,f + (AffAl ~ (Affida - A,f 
+ AAaAB + ^eij{2AiAAj + AA^AAj) 

- AfAfAo{x~AAo)+4fAf{d,x~A,)AA, . (20) 
The small corrections to the fields has to be calculated from following equations obtained by linearisation of full field equations. 

V'(^) + 4p(2 - 3p)(^) + (dk lnp)[dki^) + SkiiAAi - diAx)] - 2{p - l)(A^o - x) = J^^^P^ ^^^^ + 7^ 

V^Ax ~ dkAAk - 2{2p - 1)^ + {du \np){ekidi^ + d^Ax - AAu) - X - Ako + Ax - 2(1 - p)^ 

ekidkAAi-Af{\- p)Af + 2p{x- AA^) = -2pAx 

dnAAo + 2dnp^ + 2penkidkAx-AAk)-dnX-Snkdkij) = AAn ■ (21) 

The equations were simplified with a use of static field equations satisfied by background fields. Now the crucial observation is that 
in the Bogomol'nyi limit any forces exerted on vortices must be zero for vanishing velocities so accelerations are at least linear in 
velocities 

A;^=</(A)Af . (22) 
u for given p is a matrix in indices A, B. This relation can be iterated to give 

73 a AB 

" xA _ AB(^x-<B I ^^pq \B {C AB .BC \C /9o^ 

where once again we have preserved only terms linear in velocities. We can also make following definitions and approximations 

/ = ff^p'^p J 

/ ~ ff^p^p,l^\ J 

A/ ^ fs^X^ , 

Af « /«^Af , 

A/ « /.>,r-.TA^ 
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'pq ^'q ) 



AAk ^ At^X^ , 
AAo = a^X^ , 

With these formulas and with the gauge Ax = Eqs.(pTl) can be rewritten as 

X^lV's^ + 4p(2 - 3p)s^ + {du \np){dus^ + SuiAf^) ~ 2{p - l)(a^ - ^)] = /i^<^Af + .^.^^''Af 

X^[~d,At^ 2{2p l)h^ + {du \np){eMdis^ ~ A^^] = " < " 2(1 - P^K'^Af 



(24) 



X^[e^iduAf^ -^p{l-p).4 + 2p{^-a^)] = 



dX^ 

X^[dna^ + 2dnps^-2penkA^n-idn^+enkdkh^)] = A^Pcj^fxf . (25) 
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They should yield field deformations in approximation up to terms linear in velocities. In the limit of very small velocities field 
deformations become small as compared to background fields so our linearisations of field equations with respect to deformations 
are justified. Because full field deformations are regular functions then in the limit of slow motion also Eqs.(p5|) must have regular 
solutions being good approximations to full deformations. With this in mind we can take for granted existence of regular solutions 
and derive necessary conditions for the regularity. The unknown parameters in Eqs.(^5|) are elements of matrices ujp^{X). Our 
strategy from now on is to adjust such unique values of these parameters which allow solutions to be regular. Once we know the 
parameters we will also know equations of motion for A's linearised in velocities 

A^=c^;l,^Af . (26) 

Since we know linear part of the effective Lagrangian (^8|) the knowledge of these equations enables us to restore also its quadratic 
part up to total derivatives. The general form of the Lagrangian is 

Leff = <^(A)A^Af + b^f{X)X^X^ . (27) 

The metric tensor on the moduli space g must be symmetric under exchange of pairs of indices {A,p) and (_B, q) and invertible 

ig-Xfafr'' = • (28) 

Equations of motion linearised in velocities are 

AB"A AB BA ■ A ^^^^ ^^"^^ A'D 

'^9pq \ + {bpq -bgp )X.p + {-Qj^ ^)^)\-^s • (29) 

The components of the metric tensor must solve following equations 

nBC, ,CA . ^ (,AB ,BA\ , / '^"rg I_!Z_UC - C^n) 
9qr ^rp + ^^0^^ " V ) + 2^~dxf ^ "Q)^'^-^ ~ 

This is a system of 2n x 2n linear inhomogenous equations. The basic condition for the system to have an unique solution is that 
matrix lo is invertible in pairs of indices {A^p) and (C, r), detuj ^ 0. This means (see Eq. (^6|) ) that whatever is the small velocity it 
is always a source of acceleration already in linear terms. It should be a quite generic case except some anomalous sets of measure 
zero in a model with magnetic interaction. The metrics can be extended to these exceptional points by continuity. One of such 
points is certainly the limit of infinitely separated vortices where magnetic interaction degenerates to a purely topological term 
(Eq. ( |l8| ) . But in this limit the effective Lagrangian can be accurately calculated with a help of the product Ansatz of independently 
Lorenz-boosted vortices. The quadratic term reads 

V 

The second and third terms in Eq. (|30|) are explicitely antisymmetric under exchange of pairs {A,p) and (B, q) so the first term also 
has to be antisymmetric lo^ g = —gui. This condition means that to leading order acceleration is orthogonal to velocity with respect 
to metric g. Once again it should be so for forces due to magnetic interactions and they are the only forces linear in velocities. If 
this condition is satisfied we are left with 2n^ — n independent equations necessary to establish the same number of metric tensor's 
components. 

4 Two vortices in center of mass frame 

We will consider the by now classic example of two vortices in CM frame. The system is well described by two parameters Ai and 
A2 which can be identified with former notations 

Ai = Aq , A2 = Aq . (32) 

The Higgs field is c/) — {z^ — X)W{z, z, A). Eqs.(|25|) are on both sides linear in velocities so they can be cast in a form of a velocity- 
independent matrix multiplying a vector of velocities. For the product to be always zero any velocity vector must belong to kernels 
of the matrices. A matrix is zero if and only if it anihilates any vector from the basis spaning the space of velocities. If we choose 
Ai 7^ and A2 = Eqs.(p5|) reduce to 



V^sd) + 4p(2 - 3p)s(i) + {d, lnp)(5,.(i) + suA^^^) - 2{p - l)(a(i) ~ ^)] ^ h^^^c.^' + s^^h 

-afeA^ - 2(2p ~ l)/i(i) + (dk lnp){skAs^'^ - aW) = - a(^) - 2(1 - p)s^^^]co^' 

£fci9,^|^^-4p(l-p)s(i)+2p(|f -a(i)) = 

OAl 

5„aW+2a„ps«-2pe„fc4'^-(9„^+£„fe9fc/i(i))] = A^^^^^^ . (33) 
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On the other hand for the choice Ai = and A2 7^ we obtain 

V^sf^) + 4p(2 - 3p).(2) + (a, lnp)(9,.(2) + ^^^^(2)) _ 2(p - l)(a(2) - ^)] = /i(^)c.-42 + s(^^u;^^u;^' 

- 2(2p ~ + lnp)(efc,azs(2) - ) = - a^^) - 2(1 - p)s^^^]co^^ 

These two sets of equations have to be satisfied simuhaneously. 

Now let us consider the limit of coincident vortices. Such a configuration is rotationaUy symmetric and it should be all the same 
what is the direction in which it is splitting. This motivates the limiting form of the effective Lagrangian 

ieV -ho^\X |« S^'^XaXb - goe^'^XAXB , (35) 



with the coefficient ho and the power ^ to be determined, go can be extracted with a help of Eq.(17) and fluctuations ( |ll| , |l5| ) 



(37) 



A value of go was estimated numerically to be 50 ~ 0.0194. 
Motivated by null powers of A's in Eqs.(B3,R3h and by 



dx sin2(? dx cos 20 
dXi ' 8X2 



in the limit of vanishing A's we can restrict field deformations to following forms in polar coordinates r, 9 

s^^) = s(r) sin 26* , s'^^ = -s(r) cos 26* , 
a(i) = a(r) sin 26* , a*^' = -a(r) cos 26* , 
Al^^ = b{r) cos2e , = 6(r)sin26' , 

A^i) = c(r) sin 26* , A^^^ = -c(r) cos 261 . (38) 
together with a form of the matrix uj"^^ — Loe^^ . Eqs.(^3|) and ( ^4|) reduce to 

s" + ^-^+w2s + 4p(2-3p)s + ^(s' + c)-2(p-l)(a--^) ^ uoH^ , 

6' + - + --2(2p-l)H2 + -(&-— ) = -^+u:a-2uj{p-l)s , 
r r p r r'' 

c' + ^ + y-4p(l-p)s + 2p(-^-a) = , 



2 2H2 

a' + 2p's -2pc + ujb = — - -\ , 

2a 2 , 

— + 2pb-ujc = — + H'^. (39) 

LD is an adjustable parameter we have to choose in such a way that solutions are regular. Short inspection shows that fourth equation 
in the above set of equations can be derived from the fifth, second and third. Regularity of A/ means that s can not be more 
divergent then 0(r~^). Thus we can expand the regular solution around r = as 

00 CXD OC 00 

s(r) = ^ sur^ , a{r)=^akr^ , 6(r) = ^ 6fcr'= , c{r) = Y,Ckr'' . (40) 

k=-2 k=0 k=l k=l 



Substitution of first few terms in the expansions to Eqs.(|39D shows that they are in contradiction unless uj — —2. If we adopt this 
value of u) following leading terms will be obtained 

s = {-v-iY^..., 

a = (-u - 2^)r^ + ... , 
b = {~2v ~2^)r + ... , 

c = (2i; + 2e)r , (41) 
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where v ~ 1387000 is a coefBcient in the expansion H2 = + vr^ + ... and ^ is a free parameter coming from a "homogenous" 
part of the solution. For the solution regular at infinity numerical analysis has given ^ « —1387017.5. 

With the value of w = — 2 we can conclude that the limiting form of the effective Lagrangian for A ^ is 

L^J^f' = g,{]^5^^XA\B ~ s^^XaXb) ■ (42) 

In a generic case of two vortices in CM frame at [R, <d) and (i?, + tt) the effective Lagrangian must take a form 

Leff ^ F{R)if + 2G{R)RRe + H{R)R'^e^ + B{R)Re . (43) 

The functions F, G, H are to be determined. It is a general form of the metrics tensor invariant with respect to rotations. The 
function G is in general nonzero since the Chern-Simons term breaks parity invariance. Eq.(|2^) also contains terms which break 
parity and there does not seem to be any reason why these terms should vanish. In polar coordinates the w-matrix reads 

AB ^ -h\RB{R)] ^ ^ JV^.-l.AC^CB (44) 

with J(i?) being the total spin of two vortices separated by a distance 2R. It is an invertible matrix and acceleration is indeed to 
leading order orthogonal to velocity. I have attempted calculating uj but because there is less symmetry (less constraints) in the 
problem for generic R then for i? = it can not be extracted just from the asymptotices close to zeros of the Higgs field. Matching 
with asymtotices at infinity would be necessary. 

Finally a comment on analogous scattering of vortices in Abelian Higgs model (AHM) is in order. In this model G{R) = and also 
B{R) — 0. Similar considerations as just above lead to a conclusion that uj^^ — 0. Thus in AHM accelerations are at least quadratic 
in velocities. In equations analogous to ( pT| ) terms linear in acceleration and its time derivative can be neglected as compared to 
those linear in velocities. These terms were indeed neglected in Appendix B of Ref.[Q and it was shown that rearrangement of 
effective Lagrangian analogous to ( pO| ) with a help of equations fulfilled by deformation leads to the same expression as that derived 
by Samols Here we have shown justification of these steps when performed in AHM. 

Let us consider a general Bogomol'nyi theory and try to decide what are the conditions for the effective Lagrangian to be purely 
quadratic in time derivatives. A Lagrangian of a theory can be written as 

L = Gabm^ai'b + Ka^i'a - e[^P] (45) 

where i/i's are a set of fields, e[ip] is a static energy density functional and Ga6['0] is an invertible, symmetric and positively definite 
tensor. The only contribution to the linear part of the effective Lagrangian is 

L^^f) ^ K4i;]^a , (46) 

where like in all of this paper ip^s are the fields of the static self-dual background and the time derivative means a total derivative 
with respect to time-dependent parameters. This term certainly vanishes if ifa[V'] = for the given background. More subtle 
possibility is that the whole expression (pq ) can yield zero result when its spatial dependence is integrated out. Relativistic gauge 
theory contains linear terms like in Eq.(|45|). In the Abelian Higgs model such a term is like —diAoAi + etp^i^Aox but the self-dual 
background has the property that = and that is why there is no linear term in the effective Lagrangian. An interesting example 
of Maxwell-Higgs self-dual model with relativistic kinetics can be found in pT| . A uniform background charge density forces nonzero 
Aq. Vortices in this model feel both Magnus force and mutual magnetic interactions. 
The quadratic term of the general effective Lagrangian is 

L'f)f = Gabmi^'aA + ^^'aA^b + A^a^^b) + ^ii>aA^b " ^Mb " A^aA'^,) - .J^[^]A^^A^, , (47) 
"JJ dipb HaOiPb 

where we have introduced time derivatives of background fields and deformations of fields. Now we restrict to the case of L^^fj = 0. 
This means that accelerations are at least quadratic in velocities and the above formula can be reduced to 

= Gab^Mb + ^(V^aAV'fa - ^Mb) - y^^iV-lAV-aA^b , (48) 

We can linearise field equations with respect to velocities. The field deformations must satisfy a simple equation 

AV'fc = -T^)Wh ■ (49) 



SipaS'ipb Slpb Slpa 

This relation enables us to simplify Eq.(^8|) to the following compact form 

Leff = Gab[lp]li'ai^b (50) 

plus higher order terms negligible in adiabatic approximation. 

Thus whenever there is no linear term in the effective Lagrangian the quadratic term can be correctly calculated just with a 
help of background fields with their parameters promoted to the role of collective coordinates. 
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5 Conclusions 



The limiting form of the term quadratic in velocities of the effective Lagrangian for two Chern-Simons vortices was extracted 
from equations satisfied by deformations of the fields with respect to static background. This form shows that as we trace locally 
trajectories of vortices passing one over another there is the celebrated right-angle scattering. Globaly if two vortices were pushed 
from a large distance one against the other with zero impact parameter they would avoid direct collision their trajectories being 
curved by charge-flux interactions. A difference between a total spin of the pair of vortices when they are infinitely separated and 
when they sit on top of one another is —2tt (for k ~ 1). Thus the necessary condition for the zeros to meet and the right-angle 
scattering to occur is that an impact parameter d with respect to the center of mass and an initial velocity v satisfy dv = This 
condition can become sufficient only for d small enough because vortex-vortex magnetic interactions are falling exponentially with 
a distance. 

The local right angle scattering is a hint that moduli-space manifold is similar to a smoothed cone. The missing volume can 
show itself in a thermodynamics of a vortex gas by an excluded volume in van der Waals state equation similarly as for the vortex 
gas in Abelian Higgs model |^ . 

Finally let me stress that problems partially overcome in this paper are not at all special to relativistic Chern-Simons models. 
They can also appear in nonrelativistic and nonvariational models so celebrated in Condensed Matter Physics because Schrodinger 
or diffusion terms are linear in time derivatives on one hand and on the other hand the effective action method may be not able to 
describe the whole variety of dynamical phenomena. I would like to address such problems in a near future. 

Acknowledgement. I would like to thank Prof. Henryk Arodz and my coUegues from Department of Field Theory for discussions 
during the preparation of this work. Special thanks should be directed to Prof. Ki-Myeong Lee for his critical remarks on the first 
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